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We study a 17(1) x (7(1) system with two species of loops with mutual 7r-statistics in (2+1) dimensions. We 
are able to reformulate the model in a way that can be studied by Monte Carlo and we determine the phase 
diagram. In addition to a phase with no loops, we find two phases with only one species of loop proliferated. 
The model has a self-dual line, a segment of which separates these two phases. Everywhere on the segment, 
we find the transition to be first-order, signifying that the two loop systems behave as immiscible fluids when 
they are both trying to condense. Moving further along the self -dual line, we find a phase where both loops 
proliferate, but they are only of even strength, and therefore avoid the statistical interactions. We study another 
model which does not have this phase, and also find first-order behavior on the self-dual segment. 



Systems with statistical interactions between particles arise 
in many contexts. For example, Laughlin quasi-particles 
in the fractional Quantum Hall Effect have mutual statistics 
which depends on the particular stated As another example, 
spinon and vison excitations in Z2 fractionalized phases in 
quantum magnets have mutual tt statisticsi^ 

While gapped topological phases are well understood, the 
phase transitions between them are less explored. One reason 
is that the character of the transitions is a dynamical question, 
while in cases with statistical interactions, the system path in- 
tegral contains generally complex phases and hence has sign 
problem in direct Monte Carlo (MC) approaches. 

References [Handle considered the toric code model, which 
can be also viewed as an Ising gauge theory with Ising matter 
fields Z In a formulation containing Ising matter and Z 2 fluxes 
(visons) as elementary particles, the two have 7r-statistical in- 
teraction, and one cannot simulate large systems in these de- 
grees of freedom. On the other hand, the matter-gauge field 
formulation does not have the sign problem and was studied 
in detail in MC. 6 This system has two phases: 7 deconfined and 
confined. The confined phase includes both the Higgs regime 
(crudely viewed as condensation of matter out of the decon- 
fined phase) and confining regime (viewed as condensation of 
visons), and there is a path connecting the two regimes that 
does not cross any phase transitions. Resembling somewhat 
liquid-gas system, there is also a finite segment of first-order 
transitions separating the two regimes along the so-called self- 
dual line, where the Z2 charges and Z2 vortices have identical 
interactions. We can think of the two species as being im- 
miscible on the first-order segment. Reference focused on 
the transition near the tip of the deconfined phase along the 
self-dual line and possible scenarios how the Higgs and con- 
finement transition lines join. 

In this paper, we study a U(\) x 17(1) system with n statis- 
tical interactions, which appears in effective field theories of 
frustrated quantum antiferromagnets^^ and other areasJiia 
We consider concrete lattice realizations that can be reformu- 
lated in a sign-free manner and explore these using MC. 

Specifically, we consider a system of two loops with short- 
range interactions and 6 = tt statistics, 
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Here J\ and J2 represent conserved integer-valued currents, 
V • Ji = V • J2 = 0, residing on the links of inter-penetrating 
cubic lattices 13 . The last term is a statistical interaction where 
the cross-linking of the two loop systems is calculated with the 
help of an auxiliary "gauge field" 3,2 whose flux encodes the 
J2 currents, J2 = V x 02. Just like in the Z2 x Z2 case, we 
can reformulate the model as a special matter-gauge system 
amenable to MC study. 

Figure [T|is our main result and shows the phase diagram of 
this model. For small ti and t%, there are only small loops 
(phase in the figure). If we fix ti and increase t\, the J 2 
loops remain small, while beyond some critical t\ coupling, 
J 1 loops condense via an XY transition (phase I); we get an- 
other phase (II) if we keep t\ small and increase ti- Unlike the 
Ising case, the two phases I and II are distinct, and the follow- 
ing question arises. Suppose we are increasing t\ and along 
a self-dual line t\ =ti = t where the two species have iden- 
tical interactions, so both equally want to proliferate at some 
point. One possibility is that they behave as immiscible flu- 
ids and phase separate. Another possibility is a regime where 
both loops are present in some critical soup, which would be 
an example of an unusual phase transition. Such a question is 
of much recent interestj&i^— 

In the present study we find that in the above specific model, 
the first scenario happens and we have first-order transition 
between phases I and II. Interestingly, if we continue increas- 
ing t\ and t2, the two loops eventually condense simultane- 
ously but only in even strength for both J] and J 2 , while there 
are no large loops of odd strength for either species. By go- 
ing into this "paired- Ji" and "paired- J2" phase (labelled III 
in Fig. [U, the loops avoid the destructive interference effects 
of the statistical interaction. 

Returning to the transition line I-H, the first-order transition 
is strongest near the two ends of the segment where respec- 
tively phase I or III opens up. We also explore what happens 
when we modify the model to eliminate phase III, which we 
realize by simply restricting | J\\ < 1, | J%\ < 1, thus prevent- 
ing pairing within each species. The strength of the first-order 
character indeed decreases as we increase t, but for all such 
accessible values the transition remains first order. 

Model and Monte Carlo Method. In order to reformulate 
Eq. ([T|i in a sign-free way, we pass from J\ variables to conju- 
gate 27r-periodic phase variables by formally writing the con- 
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straint at each r: 

<5[V • Jt(r) = 0] = / d(f> r exp[-# r (V • J x )]. (2) 

J — 7T 

To be precise in our system with periodic boundary condi- 
tions, we also require total currents of J\ and J 2 to vanish. 
In this case we can write J 2 = V x a 2 and the action (Q~|) is 
independent of the gauge choice for a 2 . We enforce the zero 
current in J\ with the help of fluctuating boundary conditions 
for the <f>-& across a single cut for each direction \i = x,y,z 

&(^J^{r)5 riltQ )= / cxp[-i7 M ^ Ji^(r)8 r ^ ]. 

r J — 7T r 

(3) 



This gives the following partition function: 

z= yi ru d ^rf[d^ 

. J —7T r J -7T 1 
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S[<t>,l,a 2 } 
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where the action is given by: 

[V x a 2 (r)] 2 



S[^,j,a 2 ]= 2J 
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^Villain is the 'Villain potential', which is obtained by sum- 
ming over the Ji variables: 



CXp[-Vvillain(<Ml)] = 2J eXP 
Ji — — 00 



(6) 



In the actual Monte Carlo, we use <f> r , 7 M e(— x, tt), 
a2/j(^)eZ, and perform unrestricted Metropolis updates. One 
can show that physical properties measured in such a simula- 
tion are precisely as in the above finitely defined model. 

In this work, we monitor "internal energy per site", e = 
S/L 3 , and compute heat capacity, defined as 



C = ((e 2 ) - (e) 2 ) x Vol, 



(7) 



where Vol=L 3 is the volume of the system. To determine the 
phase diagram, we monitor loop behavior by studying "super- 
fluid stiffness", which is defined for loops of flavor a as: 
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Because of the vanishing total current, we define these at the 
smallest non-zero q\ e.g., for p xx we used q = (0, , 0) and 
q = (0, 0, We focus on flavor 2 since it is more readily 
accessible in the formulation of (O. We also monitor gauge- 
invariant "magnetization", defined as: M = e 2l * r , which 
can detect flavor 1 condensates (care is needed interpreting M 
for different boundary conditions). 

Results for Unrestricted Currents. We determined the 
phases of this model by looking at the order parameters p 2 and 
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FIG. 1. The phase diagram for the unrestricted model. Phase (0) 
contains no loops. Phase (I) contains proliferated loops in Ji and no 
loops in J2 , while in phase (II) the variables are interchanged. Phase 
(III) contains proliferated double-strength loops in both variables. 



M. p 2 is non-zero in phases (II) and (III), and M is non-zero 
in phases (I) and (III). We found the locations of the phase 
boundaries more accurately by studying p 2 ■ L crossings. We 
took data in sweeps across the phase diagram (see Fig.Q~|i, and 
defined the intersection of the p 2 ■ L curves to be the location 
of the phase transition. An example of such a sweep is shown 
in Fig. [3] The sweep is 'vertical' in the phase diagram, i.e. 
fixed t-y. The symbols on the phase diagram are the values 
of ti at which the sweeps were performed. We studied the 
fine nature of the phases by looking at clusters formed by J 2 
loops. In phase (II), the largest clusters of J 2 grow with sys- 
tem size, and have arbitrary J 2 , so we deduced that this phase 
contains condensed loops of J 2 . In phase (III), the loops that 
scale with system size have even J 2 , so we judged that the 
condensed loops have even strength. In phases (0) and (I), we 
found no large clusters, and deduced that J 2 is gapped. The 
model with t\ = is a model containing only one species of 
loop22 Our value for the position of the (0)-(II) XY transi- 
tion (t 2 ~ 0.333...) is in agreement with prior work on this 
model^. For t\ — > 00, the Villain weight © vanishes except 
for a = 2irx (int), which enforces J 2 = V x a 2 = 2x(int). 
Therefore, at large t\ the (I)-(III) transition is a transition from 
no loops of J 2 to loops of even J 2 . One expects that this tran- 
sition is XY-like, and similar to the (0)-(II) transition, but due 
to doubled J 2 , it should occur at a t 2 value four times higher. 
We observed the (I)-(III) transition to occur at t 2 ~ 1.3 for 
large t\, in agreement with this expectation. 

To determine whether the loops are immiscible, or if a crit- 
ical state is possible, we consider whether the phase transition 
(I)-(II) on the self-dual line is first- or second-order. We stud- 
ied this by looking at histograms of e. In the continuous case, 
these histograms would be singly-peaked, while in the first- 
order case we expect to see two distinct peaks. An example 
of such a histogram is shown in the inset for Fig. [2] We ob- 
served dual-peaked histograms for all points on the self-dual 
line and concluded that the transition is first-order. One way 
to quantify the strength of the transition is to look at the dis- 
tance between the two peaks in these histograms. We plotted 
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FIG. 2. Plot of the energy difference between the peaks of the his- 
tograms, as a function of ti = ti = t on the self-dual line. Higher 
values denote a stronger first-order nature. The inset shows the en- 
ergy histogram for ti = ti = 0.6. The dual-peaked shape of the 
histogram implies a first-order transition. The histograms were per- 
formed at L — 24 with 10 6 Monte Carlo sweeps. The heat capacity 
intercepts were found using data from L — 8, 10, 12, 14 and 16. 



this peak-to-peak distance in Fig. [2] which is clearly non-zero 
for all t. 

We can also determine the order of the transition by look- 
ing at the heat capacity, where the signature of first-order is 
C growing as L 3 . Phenomenologically^, at a first-order tran- 
sition the energy distribution is described by a sum of two 
normal distributions, P(e) = c + p + + C-P-, where c+ and 
c_ are the weights representing how much time the system 
spends in each state, and p + and p_ are two normal distribu- 
tions centered at e + and e_. Computing heat capacity based 
on this ansatz gives: 
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where e + — e_ is the peak-to-peak distance, and Bis a volume- 
independent constant. We plotted ^ vs. ^ and found the 
y-intercepts of these plots, which should be equal to A. If we 
assume that c + = c_ = 0.5 as a rough estimate, then the en- 
ergy gap should be equal to 2\/A The resulting estimate from 
heat capacity is plotted in Fig. [2] Both of our measurements of 
the energy gap show that the transition is strongly first-order 
at t = 0.4. The strength of the transition initially decreases 
with increasing t, before increasing again after t w 0.6. Thus, 
in this model the two condensates are immiscible. 

In addition to finding the position of the phase boundaries, 
plots of p2 ■ L like the one in Fig.[3]can be used to study the 
nature of the (I)-(III) phase transition. We can argue that at 
high t\, the phase transition is continuous. We are interested 
in whether the nature of the transition changes to first order 
before it meets up with the transition on the self-dual line. 
In a first-order transition, we expect the values of the p2 ■ L 
crossings to increase with L. We do not observe this in Fig.[3j 
or plots at other t\. We therefore suspect that the transition is 
second-order, though our data is not precise enough to rule out 
a weak first-order transition. We also obtained histograms of 
the energy and magnetization at the phase transition located 
by the p^ ■ L plots. We found no evidence of two peaks at 
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FIG. 3. piifirain) ■ L as a function of ti at ii = 1.2. Based on this 
data, we conclude that the phase transition occurs at approximately 
t-2 — 1.168, and is second order in nature. Scans like this were used 
to determine the phase diagram in Fig. [T] Here the transition is be- 
tween (I) and (III). Each data point is the result of 5 x 10 6 Monte 
Carlo sweeps. Error bars were determined by looking at the differ- 
ence between runs with different initial conditions. 




FIG. 4. The phase diagram of the 'restricted' model. The region 
with +45° lines is inaccessible in the formulation (|5j, and the cross- 
hatched region is inaccessible in this work. The transition is first 
order everywhere on the accessible portion on the self-dual line, but 
the strength decreases as t is increased. 



L = 16, 20. Since we do not know the critical point exactly, 
it is difficult to study such histograms at larger sizes, so once 
again the data indicates a second-order transition but cannot 
rule out a weak first-order one. 

Results for the Model with Restricted Currents. From 
Fig. 12 we have seen that the first-order character of the phase 
transition on the self-dual line is strongest where it meets 
phases (0) and (III). Therefore we might expect the transi- 
tion to be more weakly first-order if we could eliminate phase 
(III). This is the reasoning behind the restricted model, where 
we only allow loops with |Ji| < 1, | J2I < 1- This has been 
done in the Monte Carlo by only changing if the result- 
ing curl satisfies the restrictions, and by restricting the sum in 
Eq. © to run only over the values —1, 0, 1. 

Figure|4]shows the phase diagram for this modified system, 
which was determined using the same methods as Fig.Q] The 
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FIG. 5. Same as figureUJ but for the restricted model. 



boundaries of phase (0) are very similar to the unrestricted 
model, which is not surprising as at the (0)-(I) transition the 
proliferating loops are mostly of strength 1 . Referring to the 
Villain potential (O, we see that if the sum on the right side 
is negative, the potential is undefined. When the sum is re- 
stricted to \Ji\ < 1, this occurs for ti > 0.72135, and this 
limits the area of the phase diagram that we can study with 
this formulation. Upon using interchange symmetry, Fig. [4] 
contains a region inaccessible in this work, indicated by cross- 
hatching. 

We now investigate whether removing phase (III) has 
changed the nature of the transition on the self-dual line. We 
used energy histograms and studies of the heat capacity to 
determine the peak-to-peak distance. The heat capacity sug- 
gests a first-order transition. For t = 0.6, 0.7, the energy gaps 
were too small to be accurately determined by studying the 
histograms. The histogram for t = 0.7 is shown in the inset to 
Fig- El We cannot resolve two separate peaks, but the distribu- 
tion has a flat top, which suggests that the transition is weakly 
first-order. In order to acquire more clearly two-peaked his- 
tograms, we studied the magnetization of the system. We 
found the peaks in these to be more easily distinguished, and 
the results clearly indicate a first-order transition. 



Discussion. We studied a lattice realization of a U(l) x 
17(1) system with 7r-statistical interactions. It was helpful to 
know the location of the phase transition between I and II from 
self-duality in this non-trivial 3D Statistical Mechanics prob- 
lem. In two somewhat different models, we found first-order 
transitions on the self-dual line, which means that when both 
loops are trying to condense, they tend to phase-separate. A 
continuous transition would be an example of an NCCP 1 U(l) 
self-dual critical point 8 . We found in the restricted model that 
the first-order transition became weaker as t was increased, 
but we could not study the model for t higher than a certain 
value. If one could find a way to study the model at high t, it 
would be interesting to see if the first-order transition contin- 
ues to weaken and perhaps becomes second-order. One could 
also explore more models asking if some short-range modifi- 
cations can produce a critical loop state. There is evidence for 
a continuous transition in SU(2) spin models , 16,17 ' 25,26 but our 
system has no analog to these. 

Our study is an example of a sign-free reformulation of a 
model with statistical interactions and the power hence af- 
forded by Monte Carlo to establish the phase diagrams and 
study phase transitions. Though we determined most of the 
phase diagram in good detail, it may be useful to get a better 
understanding of how the phase transitions join at the corners 
of the (0) and (III) phases. It would also be interesting to ex- 
plore more models with statistical interactions that can have 
such reformulations. An accessible direction already in the 
present setting is to examine the model Eq. §5$ with general 
statistical angle 9. Another interesting direction is to intro- 
duce some attraction between the two loop species, to see if 
we can achieve fermionic bound states of Ji and J2 and what 
phases can be accessed in this way. 
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